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O Abstract 
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We present a mathematical formalism describing the propagation of a completely general elec- 

> 

O tromagnetic wave in a birefringent medium. Analytic formulas for the refraction and reflection 

from a plane interface are obtained. As a particular example, a Bessel beam impinging at an ar- 

bitrary angle is analyzed in detail. Some numerical results showing the formation and destruction 

c/3 of optical vortices are presented. 
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I. INTRODUCTION 



The phenomena of birefringence was aheady well known at the time when Huygens 
studied the "strange properties of the Island crystal" pQ with a view to proving the wave 
nature of light. Newton also discussed these same properties at length in his Opticks [2], 
but to prove precisely the contrary. In any case, a full mathematical description of the 
phenomena turned out to be quite cumbersome and it is only in the last few decades that 
analytic expressions were obtained for the simplest case of a plane electromagnetic wave 
in an uniaxial crystal [3]. Generalizations to more realistic waves, such as Hermite- Gauss, 
Laguerre-Gauss or Bessel beams, have been so far restricted to paraxial (or almost paraxial) 
approximations |H El E] , and to propagations parallel [7j or perpendicular [H] to the crystal 
axis. 

The aim of the present work is to obtain the most general expressions describing reflected 
and refracted (ordinary and extraordinary) waves in terms of the properties of a beam 
impinging at the plane interface of a uniaxial crystal. Analytic formulas are presented in 
as compact a form as possible. As an example of application, the propagation of a Bessel 
beam P [10] is studied; these are electromagnetic modes that propagate in vacuum with an 
invariant intensity pattern and exhibit polarization and phase optical vortices [TT] yielding 
an electromagnetic orbital angular momentum (see, e. g., Ref. [12]). However, as we will 
show in the following, the intensity pattern inside a birefringent crystal does not remain 
constant along the main axis of propagation (unless it coincides with the crystal symmetry 
axis); moreover optical vortices can be created or destroyed. Flossmann et al. [131 E] 
considered the case of a Laguerre-Gauss beam propagating inside a birefringent crystal and 
analyzed the evolution of the polarization vortices in terms of Stokes parameters. Here, we 
apply such a study to Bessel beams using our analytical expressions, and complement it 
with a study of the corresponding topological features of phase diagrams. 

Modern experiments of parametric down-conversion use crystals that are birefringent 
besides being nonlinear. These properties are particularly important for generating entan- 
gled photons with different dynamical properties determined by their source beam. Thus, 
for example, when a Laguerre-Gauss or a Bessel beam is taken as a pump, the down- 
converted light is expected to consist of entangled photons with orbital angular momentum 
|15j . However, the anisotropy of the birefringent crystals used for that purpose prevent a 
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direct identification of the expected characteristics of the ordinary and extraordinary beams 
and, consequently, of the properties of the idler and signal photons. Although these are 
effects of quantum nonlinear optics, their precise characterization must be preceded by a 
detailed description of the classical linear effects as wc present in the following. 

The organization of this article is as follows. In Section 2, a formalism describing elec- 
tromagnetic waves of any form inside and outside a birefringent crystal is presented; the 
boundary conditions at a plane interface are applied in order to obtain explicit expressions 
for the refiected and refracted waves. The resulting equations are used in the Section 3 to 
describe a Bessel beam incident at an arbitrary angle on the crystal interface (though the 
formalism is completely general, the optical axis is taken perpendicular to the interface for 
simplicity); we show that the reflected and refracted beams are given in terms of a single 
circuit integral. Lastly, we discuss numerical results for some specific parameters of the 
incident beam, following in particular the evolution of optical vortices, both for polarization 
and for phase. 



II. PROPAGATION IN A BIREFRINGENT MEDIUM. 

Consider an anisotropic medium described by a dielectric tensor Cjj or, alternatively, a 
dyad 'e such that the electric displacement is D = e~ • E and B = /xH. For a birefringent 
medium, 

el + Ae s s, 

where s is the axis of symmetry of the medium, Ae = ey — e, and e and ey are the permeabilities 
perpendicular and parallel to the symmetry axis respectively. 

The Maxwell equations in the absence of free charges and currents are 

V-B = , VxE+ — = 0, (1) 

V-D = 0, VxH-^ = 0. (2) 

at 

It is straightforward to see by direct substitution that their general solution for a birefringent 
medium is 

E° = sxV^^, = -^V(s- V*^) +/^^^ s (3) 
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and 

= lv(s ■ V^^) - e^^ s , = s X V^^ , (4) 

where and are Hertz potentials [TB] satisfying the two equations: 

- e/i^^ + V^^^ = (5) 

and 

-e||e;u^^ + V-?- V^^ = 0. (6) 

As it is well known, there are two fundamental modes: the ordinary wave with s ■ E° = 
and the extraordinary wave with s ■ = 0. Clearly, \E'*^ and \E'^ are associated to the 
ordinary and extraordinary waves respectively. 



A. Boundary conditions 

In the following, we restrict our analysis to harmonic fields of the form \E'()f:, x) = 
e~*^*^/'(x). Let the vacuum be defined as the region 2 < and consider a wave that impinges 
from z < and propagates inside the medium, in the region z > 0, with wave vectors 
k*^ = {kx,ky,k'^) and = {kx, ky, kf), for the ordinary and extraordinary components 
respectively. 

The general solution of Eqs. (|5| and ([6]) for 2; > is 

il;'-''''^\x,y,z) = ^ [ dkxdky e^^-^^.y^^^'"'- , (7) 



27r „ 

where k'^f'^^ in the integral must be taken as a function of k^ and ky^ namely k'^ = (e/icu^ 



kx ~ kyY^"^, and kf as the solution for k^ of equation (12). Hereafter, the factor e"*"^* is not 
included for simplicity. The two-dimensional Fourier transforms in the above formula are 
defined as 

V^^^'^H^x, fc,) = ^Jdx' dy' e-^'-^'~''yy' ^(^•^)(x',y',0+). (8) 

The (three-dimensional) Fourier transforms E and H of the electric and magnetic fields 
follow from Eqs. ([s]) and (|4]). Thus 

1 



^ s-k^)k^ -^xuj' s V^^(k^) , (9) 
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and 



H 



o 



- -(s • + euJ" s V'°(kO), = a; s X V^lk"") , 

/i J 



in obvious notation, with the conditions 



e/io;^ - (kO)2 = 0, 



e ey/x a;^ - k^ -e-k^ = 



(10) 

(11) 
(12) 



B. Reflection and refraction 



In order to study the reflection and refraction of the wave, we write the electric vector E 
in vacuum (that is, for 2; < 0) in the form 

1 



27r 



dk^dky e'^^'^+'^yy 



(13) 



where now 

K = {u^' -kl- klf/^ (14) 

and 'Ei''^'^\kx, ky) are the two-dimensional Fourier transforms of the electric field components 
of the incident and reflected waves, E(^'^)(x, y, 0~) at the interface; similar equations apply 
to the magnetic field component. 

The boundary conditions imply the continuity of Ex, Ey, and Hy at the interface 
z = Q (the continuity conditions on and arc not independent since, from the Maxwell 
equations, iujDz — dyH^ — dxHy and iuBz — dxEy — dyE^). It is convenient to express each 
Fourier transformed component of B and E^ in the vacuum region in terms of only E^ and 
Ey using the Maxwell equations. For the incident field: 



^- l{hEi + kyEl) 



k^u 



KkyEi + {kl + kl)El 



kzco 



(kl + kl)Ei + kxkyEl 



Bi = ^[-kyEi + kxEl). 



(15) 
(16) 
(17) 
(18) 
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These equations can be rewritten in diad notation as 

X = -k,uj (cj^l - kxkx)"^ (19) 

(here and in the following, V_l = {Vx-, Vy) for any vector V). 

For the reflected field, it is only necessary to change the sign of k^. Accordingly 

e, X (B^ + B^) = -k,uj {uj^l - k^k^)"^ (Ei - E^) , (20) 

and the boundary conditions take the form 

Ei + E^ = E2 + Ef , (21) 

and 

Ei - Ei = - {uj'l - k^k^) [e, X (H« + H^)] , (22) 

where E*^, E^, H*^ and H*^ are to be taken from ^ and (10). The above equations form 
a set of four equations for the four unknown functions E^, E^, ip^ and ip^ in terms of El. 
and Ey. Explicitly we have 

Ei + Ei = P + Q (23) 
Ei - Ei = R + S , 



where 



and 



P = cu(s X k^)x , Q = (s ■ k^') kx - ^iUJ'' (24) 



R = -^[yu (s ■ k*^) X kx — e c<j^ x s + e ■ (s x k^) k^], 

kz 

S = ^(s,k^ - fcf s^) + ^[ (s ■ k^)fcf - s.ki] k^ . (25) 



Therefore 



r^o ^ 2 [ez X (Q + S)] ■ E^ 
^ (Q + S)-[(P + R)xeJ ^ ' 



TE o X (P + R)] ■ 



^ - ^ (Q + S)-[(P + R)xe,] ' ^^^^ 
for the transmitted fields, and the full refracted electromagnetic field is given for all its 
components by Eqs. ^ and (10). 
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Also 

E^ = l(P-R)V^o + ^(Q-S)^^, (28) 
Ef = l-^[KE^ + kyEf), (29) 

for the reflected fields. 

The electromagnetic potentials ip^^'^\x, y, z) can be obtained by Fourier transforms from 
EqJ7|, from where the transmitted electric fields E'^*^''^) are obtained with Eqs. \m and 



Similarly, the reflected wave follows from the Fourier transform of Eq. (28). 



C. Crystal axis perpendicular to interface 



The above equations can be solved for a general orientation of the crystal axis, although 
the resulting expressions are somewhat cumbersome. They do simplify considerably in the 
particular case of the crystal axis perpendicular to the interface. Accordingly, if s = e^, then 

e\\^^u:^ = kl + kl+'Mkff (30) 
" e 

and therefore 

k^ 



Q±S = A;,(-^±l) kj 



(31) 



from where it follows that 



P±R = cu(l± ^ ) e, X ki 
jjk^ 



(Q + S)x (P + R)] =k,kluj(^'^ + l^ (l 



1.0 

fJ'kz 



Thus 



1 



iP''{x,y,z) = - I dk^dky e"'-''+"'yy+''' 

TT 



klifik, + kO) ^ ' 



ip'^{x,y,z) 



1 



TT 



dk^dky e'^-''+'^yy+'^^^ 



kl{eK + k^) ^ ' 



(32) 



(33) 
(34) 



and the electric field components of the ordinary and extraordinary waves can be obtained 
from Eqs. ^ and Q. Notice that we have used the Maxwell equations in order to substitute 
(e^ X k_|_) ■ E-'^ = ojBI and k_|_ ■ E^ = —k^El in the above equations. 
As for the reflected wave, 



E^ 



27r 




dkjky e^^-'^+^^yy-^^^^ E^ 



(35) 



where E^ is given by Eqs. (28) and (29) 



III. BESSEL BEAMS 



As an example, let us consider a Bessel beam that impinges at a given angle onto the 
surface of the crystal. For simplicity, we consider only the case of the crystal axis perpen- 
dicular to the interface. If uj is the frequency and C, is the axicon angle, then the incident 
beam can be written in the general form: 



E 



1 



2txuj 



1 / k ■ e \ "^—1 
- (£^cosecC + iB cot C) ( -r, ^ 1 

2 Vdk-e„ 



- -i ScosecC — iB cot C) ^ 
2 Vzk-e 



+ B{-^ 



5(k- eg -cjcosC) 5(|k| - , (36) 

where is the direction of propagation of the beam and is the standard (non normalized) 
left hand polarization vector perpendicular to e^. The magnetic field B is given by the same 
expression as above with only the interchange B ^ £ and £ —B. 
Taking = and = e^. + iey^ we obtain: 

1 



E 



-iuit+ikzZ 



2k, 



{lu£ + ikzB)Jm-iik±p)e 



+ {u£ - iKB)Jm+i{k^p)e'^'^^^^\ex - ley) 



+ e 



'^'BJm{k±p)e'"'^e, 



(37) 



where k^ = uj cos ( and kj_ = uj sin (, in complete accordance with the standard expressions 
for Bessel beams [9j (see, e. g., Eq. (2.6) of our previous paper [T7]). 
In order to consider a beam impinging at an arbitrary angle a, we set 

Bq = sin a By + cos a 



and 



Gp = + i{cosa By — sin a e^). 



without further loss of generality. 



Due to the first Dirac delta function in Eq. (36), the k^ integration can be performed 
directly setting 

kz = {uJ cos C ~ ky sin a) / cos a (38) 
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in the integrand. Thus the two-dimensional Fourier transform of the electric field at the 
interface takes the form 



Ezikxi ky 



Ijj 



- tan a (£^cosec(^ + iB cot Q \ ^ 



ilk ■ 



+ - tan a (£^cosecC — iB cot Q ^ 
where now 



If . p \ m+l / k ■ P \ 1 



i k ■ 



5(|k| -cj) 



(39) 



k ■ Bp = fcj: + %{ky sec a — uj tan a cos C)- 

The magnetic field component Bz{kx,ky) is obtained from the above expression by simply 
changing B ^ £ and £ ^ —B. 



The next step is to substitute the above expressions in Eqs. (33 ) and (34 ). In order to per- 



form the corresponding integral, the following change of variables to ellipsoidal coordinates 
(f/, V) is appropriate: 



kx = oj sin a sin^ coshf/ cosl^ 

ky = u sin a [cos C + sin ^ sinh U sin V] . 



It then follows after some straightforward algebra that 
k ■ Bp = CO sin ( 

and 



cosh U \ ( sinh f/ 

cos V \ %\ . , I sm V 



cosh f/o 



sinh t/o 



5(|k| — ijj) dkx dky = u cos a 6{U — Uq) dU dV 

where tanhf/o = cos a. 
Summing up, 

kxiy) = tusin^cosV^ 

ky(y) = OJ (cos a sin^ sinV^ + sina cosC) 
kziy) = ti;(— sinasin^sinl^ + cosacosC), 

kliV) = kliV) + ejV) and 



k^iy) = ^e^,u^-kliy), 

kf{V) = /e^,uj^--kl{V) 

V ^ii 



(40) 



(41) 



(42) 



(43) 



The final result for the ordinary and extraordinary waves can be expressed as a circuit 
integral: 

TT Jo 

sin a 



ixojkziy) 



X 



sinC 



kl{v)[^,K{v) + kO{v)] 

{—B COS V + S cos C sin V) + S cos aj (e^ x k) , 



and 



E^(r) = -i-^ r dV e*[r-k"(^)+-^l 



K{V) 



kl{V)[ek.{V) + kf{V)] 



X 



sma 



sinC 

For the refiected wave 



E5(r) = 



1 



^i[xk:,iV)+ykyiV)-zk,{V)+mV]_ 



X 



'0 

sin a 
sinC 



(— B cos V + £ cos C sin V) + £ cos ctj x k) 



+ UV) 



and 



fcf(\/)-6A:,(\/) 
A;f(V) + eA;,(V) 

1 . 



sm a 
sinC 



{£ cos V + B cos C sin F ) + B cos a 



27r 



27r 



A;f (y) - eA;2(y) r sin a 



kf{V) + ek,{V) 



sinC 



cos V + B cos C sin V) + B cos a 



(44) 



(45) 



(46) 



(47) 



Notice that the refiected beam is invariant under propagation along the main propagation 
axis. This axis makes an angle —a with the normal of the interface surface. 

For a beam impinging perpendicularly to the interface, a = 0, the above integrals can 
be analytically solved in terms of Bessel functions with a proper scaling of the parameter ^. 
The result is 



E' 



{nk, + kO) sinC 
+ Jm+i{k±p)e'^"''-'^He, - iey) I , 



[j„i_i(/c_Lp)e* 



("1-1)01 



(48) 



E^ 



^ k^ cos C ^-iuit+ikfz 



(ekz + kf) sm( 



Jm-i{k±p)e 



i{in—l)4>i 



+ iey) 



(49) 
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FIG. 1: Intensity pattern of an incident Bessel beam of order m = 2 impinging a calcite crystal 
at an angle a = vr/S as seen in a plane (a) perpendicular to the main axis of propagation and (b) 
parallel to the interface (notice the change of scale). The beam is circularly polarized, E = iB^ and 
C = vr/SG. The width of the observation window is measured in units of A = c/a; 



where kz = <jj cosC, = cua/ e/i — sin^ C, and fcf = oj^j e/i — (e/ey) sin^ C,. These expressions 
show exphcitly the polarizing effect of birefringence. 

Figure Q illustrates the intensity pattern of a Bessel wave in a plane perpendicular to the 
mean direction of propagation and in a plane parallel to the interface. The corresponding 
intensity patterns of the reflected, ordinary and extraordinary beams for a calcite crystal 
(e = 2.748964, ey = 2.208196) are illustrated in Figs.(|2jj3|. In these examples, the second 
order (m = 2) incident Bessel beam is taken at almost the paraxial limit ^ = 7r/36, the 
incidence angle is a = 7r/3 and the incident wave is a linear superposition of a transverse 
electric and a transverse magnetic beam with equal amplitudes and a 7r/2 phase difference: 
S = iB. The intensity pattern of the reflected beam (Fig. 2), at the usual angle = —a, 
is similar to that of the incident one. As mentioned above, it is invariant under propagation 
and exhibits an elliptic-like symmetry. The ordinary and extraordinary waves are illustrated 
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FIG. 2: Intensity pattern of the reflected beam resulting from an incident Bessel beam of order 
m = 2 with the same configuration as described in Fig. 1. The plane of observation is tilted to be 
perpendicular to the main direction of propagation. This beam has an intensity pattern invariant 
under propagation along that axis. The width of the observation window is measured in units of 
A = c/uj. 

in Figs. (3); for both refracted waves, it is possible to identify an axis of symmetry, though 
the beams are not propagation invariant. A textbook result is that for a plane wave, the 
refraction angles of the ordinary and extraordinary waves, and respectively, are 
given in our case by 

• PW sin a 

sma^ = (50) 



• PW sma 
sma^; = , . (51) 

'e/x + (Ae/ej|) sin^ a 



However, for a Bessel beam, the incident wave is a superposition of plane waves propagating 
in a cone with an aperture given by the axicon angle (; it is only in the limit C ^ that the 
above equations define the correct angle for the axis of symmetry ao and a^. In general, 
no analytical expression for an average value of ao and aE could be found. 

The intensity patterns (la) and (2a) in Fig. (3) are given at the interface surface, and all 
the other patterns are evaluated on planes perpendicular to the main propagation axis of each 
wave. The latter were obtained performing the corresponding rotations for the observation 
points and electric fields. The intensity patterns of the ordinary and extraordinary beams 
have a structure similar to that of the reflected wave near the interface, but this structure 
gets gradually more complex as they propagate. 

Since in this example the ordinary, extraordinary and reflected electric fields are ap- 
proximately contained in a plane perpendicular to their main direction of propagation, the 
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(la) (lb) (Ic) (Id) 




FIG. 3: Intensity pattern of the ordinary (1st row) and extraordinary (2nd row) beams resulting 
from an incident BcsscI beam of order m = 2 y^^th the same configuration as described in Fig. 1. 
The beam intensity pattern shows a richer topological structure as it propagates. The width of 



topological structure of their polarization can be studied using Stokes parameters [H]; these 
are defined as 



So — l-E'xP + \Ey\'^ 

C _ I p |2 I 771 |2 

5*2 = E^Ey + E^E* 

Ss = -t{E:Ey - e^e;) , 

(52) 

and are given in terms of light intensities for different orientations of an analyzer. The condi- 
tion 5*3 (x) = corresponds to local linear polarization at the point x, and 5'i(x) = 5*2 (x) = 
to circular polarization. A polarization singularity corresponds to 5*1 (x) = 5*2 (x) = 6*3 (x) = 
0. Since Sq = Sf + S2 + 5*1, this condition gives a zero intensity point where polarization is 
not well defined. 

Performing such an analysis on the incident, refiected and refracted beams, we found ex- 
tended regions where the refiected and refracted beams are approximately linearly polarized 
if the incident beam is circularly polarized in the sense that S ^ ±iB; this is illustrated in 
Figs. (4-6). As expected, the polarization of the ordinary beam is orthogonal to that of the 
extraordinary beam. The behavior of the polarization singularities is quite interesting: for 
the incident beam, an optical vortex is located at its center so that E(0) x eq = 0. But 
neither the refracted nor the refiected beams are null at that point: this means that a po- 
larization vortex may not survive the refiection and refraction process. However, 6*0 = at 
other locations both at the interface and along the produced beams and, in fact, the number 
of these singularities increases as the refracted beams evolve. This phenomenon is similar 
to that described in Ref. |T31 [E] for a Laguerre-Gaussian beam. 

Given the analytic expressions, the phase patterns can be easily analyzed for each compo- 
nent of the electric field. Such an analysis is essential for the understanding of the mechanical 



properties of the beams. For a Bessel beam of order m, Eq. (37) implies that the compo- 
nents Ei are superpositions of fields with phase m'0 -|- kzZ, where m' = m, m ± 1. The linear 
momentum density along the z direction is proportional to and the term proportional to 
determines the orbital angular momentum density along that axis, which is given by the 
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FIG. 4: (Color on line) Intensity, polarization and phase patterns for a Bessel beam of order m = 2 
with B = i£ and C = 7r/36. The patterns are produced in a plane perpendicular to the main axis 
of propagation (upper row) and at an angle a = vr/S with that axis (lower row). The second and 
third columns show the density plot of the phase for the electric field along two perpendicular 
vectors on the corresponding plane. The brighter (darker) regions correspond to phases close to tt 
(— tt). The unit of length is A = c/oj. 

formula 

i=x,y,z 

Accordingly, phase vortices with a topological charge m contribute to the orbital angular 
momentum with a factor proportional both to m and to the moduli of the corresponding 
electric field components. For reference, the phase patterns of a m = 2 Bessel beam for two 
perpendicular components of the electric field in two distinct planes are shown in Fig. (|4]): 
one perpendicular to the main axis of propagation (first row) and another at an angle 
a = vr/S with the former (second row). In the first row, the standard m' = 1,2,3 optical 
charge vortices are observed, but in the second row the phase patterns exhibit an extremely 
complex topological structure (dislocations, vortices and bifurcations) on both short and 
long length scales of order l/kz and l/k± respectively. It is also worth noticing that the 
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C) 



■100 



100 



FIG. 5: (Color on line) Intensity and polarization patterns for the reflected, ordinary and extraor- 
dinary beams at the interface resulting from an incident Bessel beam of order m = 2 with the 
same configuration as described in Fig. 1. Notice the different scales. The observation window 
dimensions are given in units of A = c/cj. 

electric field perpendicular to the tilted plane is not negligible. 

The reflected, ordinary and extraordinary beams just at the interface posses, as expected, 
a phase structure qualitatively similar to that shown in Fig. 4 in the (lib) and (lie) plots. 
In particular, the reflected wave exhibits an elliptic-like symmetry with phase vortices in the 
plane perpendicular to its axis; this is shown in the first row of Fig. (§. 

Due to the anisotropy inside the crystal, the only component of the angular momentum 
that can be conserved is the one along the crystal axis [IH]; its density is given by 

^' = iT. (54) 

i=x,y,z 

where D is the electric displacement vector. If the ordinary and extraordinary beams do 
not propagate along the crystal axis, their angular momentum along their propagation axis 
is not conserved. This in turn should be manifested as the creation and annihilation of 
phase vortices as the beams propagate in the crystal. In the second and third rows of 
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FIG. 6: (Color on line) Intensity, polarization and phase patterns for the reflected (1st row), 
ordinary (2nd row) and extraordinary (3rd row) beams resulting from an incident Bessel beam 
of order m = 2 with the same configm'ation as described in Fig. 1. The plane of observation is 
perpendicular to the main axis of propagation of each beam and is located at a distance 10500A 
from the interface plane. The brighter (darker) regions correspond to phases close to vr (— vr). Unit 
of length is A = c/iv. Notice the different coordinate scales. 

Fig. ([6]) the phase diagrams for the ordinary and extraordinary beams are illustrated at 
planes located deep inside the crystal: both beams have a rich topological structure but 
the characteristic length is much smaller for the extraordinary wave. Notice that for both 
ordinary and extraordinary waves, there is a strong correlation between the phase diagrams 
and the intensity patterns. 

IV. CONCLUSIONS 

In this work, explicit formulas have been obtained that permit an analytic treatment of 
birefringence in uniaxial crystals. These expressions are valid for arbitrary orientations of 
the incident beam and the axis of the crystal with respect to the plane defining the interface. 
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An illustrative application has been worked out in detail for a Bessel beam and a particular 
orientation of the crystal axis: it has shown that the electric fields of ordinary, extraordinary 
and reflected waves take a relatively simple analytic form in terms of a circuit integral. 
With these results, it was shown that the polarization and phase vortices of the ordinary 
and extraordinary beams evolve into a complex structure for sufficiently thick crystals. 
This structure is quite different for the ordinary and extraordinary beams and necessarily 
leads to different mechanical properties of the field. In was also shown that the refiected 
beam has elliptic symmetry, exhibits phase and polarization vortices and is invariant under 
propagation. We expect these results to be useful for the correct characterization of general 
waves in birefringent crystal as they are widely used at present in optical experiments. 
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